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Abstract 

We consider a heterogeneous cellular network with densely underlaid small cell access points (SAPs). 
Wireless backhaul provides the data connection from the core network to SAPs. To serve as many SAPs 
and their corresponding users as possible with guaranteed data rates, admission control of SAPs needs to 
be performed in wireless backhaul. Such a problem involves joint design of transmit beamformers, power 
control, and selection of SAPs. In order to tackle such a difficult problem, we apply t\ -relaxation and 
propose an iterative algorithm for the l\ -relaxed problem. The selection of SAPs is made based on the 
outputs of the iterative algorithm. This algorithm is fast and enjoys low complexity for small-to-medium 
sized systems. However, its solution depends on the actual channel state information, and resuming the 
algorithm for each new channel realization may be unrealistic for large systems. Therefore, we make use of 
random matrix theory and also propose an iterative algorithm for large systems. Such a large system iterative 
algorithm produces asymptotically optimum solution for the t \-relaxed problem, which only requires large- 
scale channel coefficients irrespective of the actual channel realization. Near optimum results are achieved 
by our proposed algorithms in simulations. 
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I. Introduction 

Fifth Generation (5G) mobile communication networks are expected to provide ubiquitous ultra- 
high data rate services and seamless user experience across the whole system [U), 0, which makes 
it necessary to offload huge volumes of data and large numbers of users from macrocells to small 
cells. In 5G networks, hyper-dense deployment of small cells will be a key factor to achieving better 
spatial resource reuse and tremendous capacity enhancement to macrocells. The services of small 
cells are provided by small cell access points (SAPs), which are low power nodes only serving 
local-area users. In order to obtain the small cell user data, the SAPs are connected to the core 
network via backhaul. 

As large numbers of small cells are underlaid with macrocells and more users are diverted 
from macrocells to small cells, providing fast and reliable backhaul connection between the core 
network and SAPs becomes a critical issue for such multi-tier heterogeneous networks 0, [f4ll . 
Wired backhaul, which uses copper or fiber cables, can provide high-rate data links between fixed 
stations. However, the cost to provide wired backhaul to all SAPs may be prohibitive when the 
number of SAPs is large. Moreover, certain locations that are difficult to be reached by wired access 
may restrict the universal deployment of wired backhaul. Wireless backhaul, which can overcome 
many of the drawbacks of wired backhaul, offers a cost-effective alternative lEJ—[EJ. Combined 
with renewable energy sources, SAPs with wireless backhaul can be established in a self-sustained 
“drop-and-play” fashion, which is especially important in countries lacking reliable and ubiquitous 
power supply [|T|. Compared to wired backhaul, the management of wireless backhaul resources, 
e.g., power and spectrum, is more complicated due to finite power and radio spectrum constraints. 

Recent works have proposed analysis and design methods for backhaul technologies from many 
aspects. A linear programming framework for determining optimum routing and scheduling of 
data flows in wireless mesh backhaul networks was proposed in |9]. Zhao el al. [ITOl considered 
the problem of minimizing backhaul user data transfer in multicell joint processing networks, 
where algorithms involving joint design of transmit beamformers and user data allocation at base 
stations (BSs) were proposed to efficiently reduce the backhaul user data transfer. Zhou et al. IfTTTl 
presented an information-theoretical study of an uplink multicell joint processing network in which 
the BSs are connected to a centralized processing server via rate-limited digital backhaul links 
employing the compress-and-forward technique. A similar scenario for the downlink of a cloud 
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radio access network (Cloud-RAN) was investigated in IIT2TI . where multivariate compression of 
different BSs’ signals was exploited to combat additive quantization noise. The spectral efficiency 
and energy efficiency tradeoff in a homogeneous cellular network was investigated in lfl3ll . where 
the backhaul power consumption was taken into consideration. Considering the overall network 
power consumption including the backhaul, Shi et al. [jT4| proposed schemes to improve the 
energy efficiency in Cloud-RAN cellular networks. For the Cloud-RAN backhauling, not only 
the routing of data but also the additional cost of baseband processing in the cloud infrastructure 
has been investigated in lfl2l . Ifl4ll . Wireless backhaul technologies have been discussed in |[5ll-lf8ll. 
Hur et al. [H proposed a beam alignment technique for millimeter wave wireless backhaul and 
investigated the tradeoff between array size and wind-induced movement. Lee et al. []6l provided 
several admission control schemes for multihop wireless backhaul network under rate and delay 
requirements. Flexible high-capacity hybrid wireless and optical mobile backhauling for small cells 
was investigated in 0- The energy efficiencies of wireless backhaul networks for different system 
architectures and frequency bands were compared in |[8). 

A vital task of wireless communications is to design schemes to meet the quality-of-service (QoS) 
requirements subject to given amount of resources. Resource management of wireless systems with 
QoS requirements has been discussed in IIT5l - lfT9| . A decentralized method to minimize the sum 
transmit power of BSs under given QoS requirements was proposed in lfT5l for a multicell network, 
relying on limited backhaul information exchange between BSs. Iterative algorithms to maximize 
the minimum QoS measure of users in multicell joint processing networks under per-BS power 
constraints were proposed in lfl6l . When it is not possible to meet the QoS requirements for all 
wireless stations with the given resources, only a subset of transmission links can be selected to 
be active. Zhai et al. IfTTl investigated the link activation problem in cognitive radio networks with 
single-antenna primary and secondary BSs and users. A price-driven spectrum access algorithm was 
proposed and the energy-infeasibility tradeoff was analyzed. The uplink user admission control and 
user clustering under given QoS and power constraints were considered in lfT8l . where algorithms 
relying on the (j-norm relaxation were proposed. Transmission schemes using semidefinite relax¬ 
ation and Gaussian randomization to select active antenna ports were proposed in Ifl9l to maximize 
the minimum user rate in multicell distributed antenna systems. 

In order to achieve enormous enhancement in spectral efficiency, the technique of “massive 
multiple-input multiple-output (MIMO)” is envisioned to be an important ingredient in 5G commu- 
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nication systems E01, nm In massive MIMO, the number of antennas equipped at BSs is much 
larger than that of active users in the same time-frequency channel. Hoydis et al. Il22il analyzed 
the number of required antennas in massive MIMO systems using different linear beamformers. 
They showed that more sophisticated beamformers may reduce the number of required antennas 
in massive MIMO systems to achieve the same performance. Fernandes et al. Il23l provided the 
asymptotic performance analysis of both the downlink and the uplink for a cellular network as the 
number of BS antennas tends to infinity. Huang et al. |[24l studied joint beamformer design and 
power allocation in multicell massive MIMO networks, where efficient algorithms to maximize the 
minimum weighted QoS measure of all the users were proposed. 

In this paper, we consider a heterogeneous cellular network with densely underlaid small cells. 
The core network are connected to the SAPs via wireless backhaul links, which provide the channel 
to transfer the small cell user data. A wireless backhaul hub (WBH) with multiple antennas is 
deployed to transmit wireless backhaul signals to the SAPs. Minimum data rate requirements must 
be satisfied when SAPs receive data from the WBH so that small cell users can be served at the 
required data rates. In order to improve the utilization efficiency of wireless backhaul supported 
SAPs and let more SAPs and their corresponding small cell users be served, it is desirable for the 
WBH to support as many SAPs as possible. Given QoS constraints at SAPs and power constraint 
at the WBH, we propose wireless backhaul transmission schemes aiming to admit the maximum 
number of SAPs into the network. Such transmission schemes involve joint design of beamformers 
and power control, as well as selection of subset of SAPs to be supported. We tackle this difficult 
problem by applying ly .norm relaxation to the original non-convex non-smooth problem and utilize 
uplink-downlink duality to transform the transmit beamforming problem into an equivalent receive 
beamforming problem. Based on the optimality conditions, we propose an iterative algorithm that 
jointly update the values of the primal and dual variables. Such an algorithm is fast and we prove that 
it converges locally to the optimum solution. Based on the solution of the fj-relaxed problem, we 
then iteratively remove the SAP that corresponds to the largest QoS gap until all the remaining SAPs 
can be supported. Furthermore, we provide a large system analysis of the above SAP admission 
control problem. As the system dimensions become large, we show that certain system parameters 
may be approximated as deterministic quantities irrespective of the actual channel realization. 
Random matrix theory is leveraged to transform our proposed finite system iterative algorithm 
for large systems so that only large-scale statistical information of the channel is required. As long 
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Fig. 1. A heterogeneous cellular network with wireless backhaul. 


as the large-scale channel coefficients and the QoS requirements remain unchanged, the selected 
SAPs will satisfy the QoS and power constraints almost surely. Simulations are carried out to 
verify the proposed algorithms for finite and large systems. The proposed algorithms demonstrate 
fast convergence and low computational complexity. For the finite system iterative algorithm, the 
average number of admitted SAPs is very close to the optimum results. For the large system 
iterative algorithm, the results of the selected SAPs accurately match those obtained by Monte 
Carlo simulations. 

The remainder of the paper is organized as follows. In Section [Til we introduce the system model 
and formulate the SAP admission control problem. The iterative algorithm for finite systems and its 
convergence property are presented in Section HEO The large system analysis and its corresponding 
iterative algorithm are presented in Section HV] Section |V] shows the convergence behavior and the 
simulation results in a cellular scenario. Finally, conclusions are drawn in Section EH 

Notation: we use bold uppercase letters to denote matrices and bold lowercase letters to denote 
vectors. CA/"(0, a 2 ) denotes a circularly symmetric complex normal zero mean random variable with 
variance a 2 . (-) T and (■) n stand for the transpose and conjugate transpose, respectively. C M and 
M m denote M-dimensional complex vectors and nonnegative real vectors, respectively. E {•} stands 
for the mathematical expectation. {u,} denotes the set made of u,, Vi ||x||, ||x|| 0 and ||xj | l stand for 
the Euclidean norm, the f'o-norm and the fj -norm of the vector x, respectively. denotes almost 
sure convergence. p(A) stands for the spectral radius of matrix A. a o b denotes the Hadamard 
product of a and b. For a vector q, q” 1 stands for the element-wise inverse of q. 
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II. System Model and Problem Formulation 


We consider a heterogeneous cellular network with N underlaid single antenna SAPs as in Fig. |T] 
The SAPs obtain their small cell user data from the core network via wireless backhaul. A WBH 
with M antennas is responsible for transmitting wireless backhaul signals to the SAPs. The wireless 
backhaul spectrum is out-of-band, which does not interfere with users in the network. In order to 
meet the data rate requirements for serving small cell users, the received data at the SAPs must 
satisfy certain minimum rate requirements. In a heterogeneous cellular network, users that cannot 
be served by SAPs will have to be served by the macrocell base station (MBS) directly. In order 
to improve the utilization efficiency of wireless backhaul supported SAPs and let more users be 
served by their corresponding small cells, it is desirable for the WBH to support as many SAPs as 
possible. Under given QoS requirements at SAPs and transmit power constraint at WBH, we propose 
schemes aiming to admit the maximum subset of SAPs that can be simultaneously supported by 
the wireless backhaul. 

We denote the wireless backhaul channel from the WBH to the ith SAP as hf, where h, e C M , 
Vi = 1, ■ ■ ■ , N. Linear transmit processing is applied at the WBH to deliver user data to the SAPs 
using wireless backhaul links. The transmit beamformer at the WBH for the ith SAP is denoted as 
u, £ C A! , such that 11u t 11 = 1, Vi. We denote the normalized data symbol for the ith SAP’s users 
that is transmitted via wireless backhaul as d t , where E {| d, \ 2 } = 1. The received signal at the ith 
SAP can be expressed as 


Vi = hf iu 



di 



dj + r/i 


( 1 ) 


N 

, h,^u 

M ^ 1 J V M 

j= ij'V* 

where rj l is the additive white Gaussian noise (AWGN) at the ith SAP, such that r] t ~ CA/"(0, rn) 
and rii denotes the noise variance. The transmit power for data of the ith SAP is Pi/M. We denote 
the power constraint for backhaul transmission at the WBH as P. The backhaul transmit power 
must satisfy '}2n =x WiPi/M < P, where Wi > 0 is the weight for the ith transmit power. The signal 
to interference-plus-noise ratio (SINR) at the ith SAP is given by 


SINFr = 


% Ih^u,- 


M 


+ rii 


( 2 ) 


Z^j=l,j^i M I i 3 

In order to meet the QoS requirements of its small cell users, the receive SINR at each SAP must 
satisfy minimum SINR requirement. The SINR requirements at SAPs can be determined from the 
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QoS requirements of their small cell users and they can be easily fed back to the core network^. 
We denote the SINR requirement for the Ah SAP as 7 Vi. 

In heterogeneous networks, small cells are deployed within the macrocell to serve users so that 
the number of users needed to be served by the MBS can be reduced to minimum. When the SAPs 
are densely deployed, the WBH may not be able to support all the SAPs simultaneously for given 
SINR and power constraints. The users that cannot be served by SAPs need to be served by the 
MBS. In order to let as many SAPs and their corresponding small cell users be served as possible, it 
is desirable to select the subset of SAPs with the maximum cardinality that can be simultaneously 
supported by the wireless backhaul. Since it is cheaper to build wireless backhaul compared to 
wired backhaul, we also minimize the total cost of building backhaul for small cells of the network 
in this way. Such a problem of SAP admission control can be formulated as 


mm x L 

p,{ui}.x 



(3) 


X > 0 


/ j i f j 1 f j 1 

where w = [uq,--- ,wn] , P = [pi,-- - ,Pn\ , and x = [x\,--- ,xn] G M+. Here Xi > 0 


represents the SINR gap of the ith SAP to satisfy its SINR requirement. If x t — 0 in the solution, 
it shows that the SINR requirement 7 , can be satisfied by the / 1 h SAP. The objective function ||x || 0 
denotes the lo-noxm of x, which is the number of non-zero elements in x. 


III. Analysis and Algorithm Design for Finite Systems 


The problem ([3]) is combinatorial and NP-hard due to the non-convex f 0 -norm in the objective 
function ll25ll . Approximate solutions of non-convex optimization problems can be obtained by 
applying convex relaxation iflOll . Ifl8l and replacing the f? 0 -norm in © with its convex envelop, i.e., 
the fj-norm. The fj-relaxed problem can then be expressed as 


mm ||x|| 1 

p,{u;},x 

s.t. constraints of ©. 


(4) 


Denote the users of the ith small cell as J. Assume the data rate requirement of the ith user to be rx, Vi £ J. Then, the SINR 
requirement for the ith SAP must be chosen such that log 2 (l + 7 ;) > r *- 



The problem © is still difficult to solve due to the need to jointly optimize the transmit power p, 
transmit beamformers {u,}, and the SINR gap x. The following lemma shows that the downlink 
transmit optimization problem © can be converted to an uplink receive optimization problem. 

Lemma 1: The downlink transmit optimization problem © is equivalent to the following dual 
uplink receive optimization problem 


min Xi 

q,{Ui},x 

SINRy > 


s.t. 


li 


1 + Xi 


Vz = 1, - - - ,N, 


(5) 


W* < p. 


M 

x > 0 

where n = [m, • • • , njv] T and q = [qi, - ■ ■ , qN] T € 


. Here we define 


SINR, U = 


M I ? 1 


2-^3 = 1 M 

Furthermore, the optimum solution of {u,} and x in © are equal to those in ©. The optimum 
solution of p in © has one-to-one correspondence to the optimum solution of q in ©. The 
problem © is a receive optimization problem in the dual uplink with the same SINR constraints, 
where h, is the dual uplink channel from the /th SAP to the WBH and w. L becomes the uplink 
noise variance at the ith SAP. Here qi/M represents the dual uplink transmit power of the /th SAP, 
and nf is the uplink receive beamformer for the /th SAP transmission. 

Proof: See Appendix [A] □ 

We denote the optimum solution of © as p*, {u*}, x*. According to Lemma |H the optimum 
solution of © can be denoted as q*, {u*}, x*. The following lemma shows that the optimum 
beamformer {u*} can be determined from the optimum uplink power q*. 

Lemma 2: The optimum receive beamformer {u*} of © is the minimum mean square error 
(MMSE) receiver, which can be obtained from the optimum uplink power q* as 

1 


% lu^hjj 2 + Wi 


( 6 ) 


u. 


u 


MMSE 


(q*) = 



—h h H + w I 

M 3 3 1 


Vi = 1, • • • , N 


(7) 


where £ is a normalization factor such that ||u*|| = 1. The corresponding uplink SINR in © is 


SINR“(q\K}) = | F h? 


N 

E 


—hh H + wl 

M 3 i + 1 


hi, Vi. 


Proof: See Appendix [B] 


( 8 ) 

□ 




For notational brevity, we define an equivalent channel matrix G, where 

Gij = |hfuj| 2 , Vi,j = l (9) 

and we omit the dependence of G tJ on u ; . After a change of variables q t = log(Jj), Vi, the 
problem © can be equivalently expressed as 

min x { 

q,{u *},x 

I* (ejLi s# G Ji^ j + w i) 

s.t. log- - -<log(l + Xi), V* = l ,-",N (10) 

Ei =1 n i e ^ 

X > 0. 

For fixed { iq }, the problem (ITOb is a geometric programming (GP) problem with variables of q and 
x. The standard way of solving a GP problem is using interior point algorithms, which is employed 
in software like cvx |[26l . Considering {u,;} also as optimization variables, one method to obtain 
solutions of (fTOb is alternately optimizing {q, x} and {u*} by solving (fTOl) with fixed {u*} and 
updating the MMSE receiver {u“ MSE } using the obtained q, respectively. However, such alternate 
optimization needs to solve (fTOl) using standard convex optimization software ll26l in each iteration, 
which makes it relatively slow in practice. In the following, we provide a low complexity iterative 
algorithm. 

We associate the zth SINR constraint in (fTOl) with the Lagrange dual variable u t , the power 
constraint with /i, and the nonnegativity constraint of x t with a,. The GP problem (fTOl) satisfies 
Slater’s condition. Hence, the Karush-Kuhn-Tucker (KKT) conditions are necessary and sufficient 
for the optimality of (fTOl) . The following lemma provides the optimality conditions for (fTOl) . which 
is key to our iterative algorithm that solves ©. 
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Lemma 3: The optimum primal and dual solutions of (fTOl) satisfy the following conditions 


4 

= max {is* — 1, 0), Vi = 1, • • • 

,N- 

(11) 

Mis* 

Q* 

\ ^ MGijyjiSj 

= 2^ f-i . *\r» * +/i Ui ’ 

(i+ x i )G is q j 

Vi; 

(12) 

x*)G*qt 

My/i 

N * 

= 2. M +Wi ' Vs; 

3=1, 


(13) 

n * 




n r ( h 

f-' M 

= P; 


(14) 

V* 

> 0; is* > 0, Vi 


(15) 


The optimum downlink power p* of © can be obtained directly from the optimum primal and 
dual solutions of (fTOl) as 


= M^is* 

M (1 + xi)G ii q!fj,*' 


(16) 


Furthermore, we have 


(1 + X*)GuP* _ GjjP* 

M 7i “.f-C. M 

3=1, 

N * 

= p 

^ M 


(17) 

(18) 


i =1 

Proof: See Appendix O □ 

Corollary 1: The optimum downlink transmit power p* of © and the optimum primal and dual 
solutions of m satisfy 


N 


a* = 


E ^+n i 


\3=1,3& 
N 


M 

Myiis*Wi 


V*Qi 

M 


, Vi = 1, ■ ■ ■ , AT; 


^ (1 + x*)G ii q k i P' 


a + , n =SINR - (q *' K}) = >■ 


N a* 

X M hjh ? + I h - V/ - 


(19) 

( 20 ) 

( 21 ) 


Here (fl9l) is obtained by substituting (fT 6 l) into (fl2l) . Eq. (l20l) is obtained by substituting (fT 6 l) into 
(fl 8 l) . Eq. (ITTI) is obtained by (fl3l) and Lemma [2] 

Based on Lemma [3] and Corollary [Q we propose Algorithm [Q that iteratively updates the values 
of the primal and dual variables to obtain the optimum p*, {u*} and x* of ©. In Algorithm Q] 
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Algorithm 1 Finite system iterative algorithm to solve (0J) 


l: Initialization: > 1 and g* > 0, Vi, such that YliLi = MP. 

2 : repeat 

3 : Store old values of q 

q = q. 


( 22 ) 


4 : Update 


Qi = 


M'-fi 


max(z/j, 1) ■ ( hf + w ^) h * 


Vi = N. (23) 


5 : Normalize 


MP_ 

q = -^q- 

n J q 


(24) 


6 : Set 


q= 2 v ^)- 


(25) 


7 : Calculate u, 


U, = 


7 ( E Jj h > h " + I h - Vz - 

\7=1 ,3 Vi 


(26) 


8 : Calculate the equivalent channel 

Gij = | hf Uj | 2 , Vi,j = !,-■■ ,N. 


(27) 


9 : Calculate /i 




max(i/j, 1) • GnQiP' 


(28) 


10 : Calculate the downlink power 

_Pi_ _ M'jjUj 

M max(i/j, 1) • Gagin' 


li: Update z/ ?: 


U = 


E 


\j=hi ^ 


GjjPj 

M 



m 

M’ 


Vi. 


(29) 


(30) 
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we first store the old values of q in (l22l) . The fixed-point iteration (l23l) and normalization (1241) are 
obtained according to (1211) and (fl4l) . respectively. Here (1251) is to ensure contraction mapping for 
the algorithm, which will be made clear in the proof of Theorem Q] After obtaining the uplink 
power q, we calculate the corresponding MMSE receiver and the equivalent channel in (l26l) and 
(l27l) according to © and ©. respectively. The value of /i in (l28l) is obtained according to (l20l) . The 
corresponding downlink power ([29b is calculated according to ([T6l) after getting the value of /i, and 
the value of z/j in ([30b is updated according to (IT9b . In this algorithm, the optimum solution of p* 
is obtained directly in (|29|) . There is no need to perform the uplink-downlink power mapping |[27ll . 
The convergence property of Algorithm [I] is shown in the following theorem. 

Theorem 1: Starting from an initial point that is sufficiently close to the optimum solution of 
©, Algorithm Q] converges to the optimum solution of © that satisfies the KKT conditions. 

Proof: See Appendix [D| □ 

Theorem Q] shows that Algorithm Q] converges locally to the optimum solution of ©. However, 
its range of convergence cannot be obtained from Theorem |T] and the result may depend on the 
initialization point. Whether this algorithm can converge globally is still an open problem and is 
left for future work. In our simulations, we observe that its range of convergence is large and even 
random initialization will converge to correct results. 

A. Connections With Max-Min SINR 

If all the SAPs can be supported, x* = 0 in the solution of (ITQb . The values of //,, Vi, decrease 
monotonically towards zero in Algorithm [Q In this case, the updates of power in (l23[)- ([24l) become 
the power iteration steps in the max-min SINR algorithm of ll28ll . Algorithm Q] will still converge. 
The output p and {u,} of Algorithm Q] are the power allocation and beamformers that maximize 
the minimum SINR of all SAPs in the system. 

B. Iterative SAP Removal 

Due to the convex relaxation, the solution of © may not be always optimal for the lo-noxm 
optimization problem ©. Therefore, © cannot be simply used as the substitution of © and we 
need to refine the selection of users based on the solution of ©. Since the solution of Xi in © 
represents the gap of the ith SAP’s SINR to satisfy its SINR requirement, it is natural for us to 
select those SAPs with small x, t . We propose to iteratively remove SAPs with decreasing values 
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of x, until the remaining set of SAPs becomes feasible to satisfy the SINR and power constraints, 
i.e., until ay = 0, Vz. In this way, we obtain the final results for SAP admission control. 

C. Admission Control Between SAPs and Users 

The data rate requirements at the admitted SAPs can be guaranteed by the wireless backhaul 
with the proposed SAP admission control algorithm. After the admitted SAPs received data from 
the core network, they need to transmit those data to their corresponding small cell users. Existing 
coordinated multipoint (CoMP) transmission schemes can be applied Ifl5l . [|T6il . However, due to the 
inter-SAP interference, there exists possibility that the users of the admitted small cells cannot be 
simultaneously supported with their given SINR requirements. If this happens, we need to perform 
user admission control within the admitted small cells in order to let as many users be served by 
their corresponding SAPs as possible. We briefly discuss the user admission control as follows. For 
simplicity of discussion, we assume each SAP serves one user. Multiple users per SAP case can 
be extended in a straightforward manner. 

We denote the set of admitted small cells as S and the power constraint for the /th admitted 
SAP as Pi, where i G S. The channel gain between the jth SAP and the /'th user is denoted as g tJ , 
where i,j G S. Within the set of admitted small cells, the user admission control problem can be 
formulated similar to ©. After ^-relaxation, we need to solve the following GP problem 

min JT x i 

p,x 

QiiPi N 7* w ■ q 

s.t. ^-■-> —-, Vz G S 

^jeSj&9ijPj + rn l + Xi ( 31 ) 

Pi < Pi, Vz G S 

x > 0, 

where Pi and ay denote the transmit power for the /th SAP and the SINR gap for the /th user, 
respectively. The problem (ITH) is similar to (flOl) . with the only difference in the power constraints. 
The problem (|3TI) with per-SAP power constraints can be solved by solving a series of weighted 
sum power constrained problems employing Algorithm Q] Following similar discussion as Il29l . Ii30l . 
it can be shown that the optimum values of ([Til) is equal to the optimum values of the problem 
max{ Wj > 0 } / ({ «;,}), where /({'«;,}) with { «;,} as the parameter denotes the optimum objective value 
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of the following weighted sum power constrained problem 

f({wi}) = min 

p,x 

S.t. 


where w l > 0, Vi e S, are the weights of power and also the parameters of /({«',;}). The 
problem (l32l) can be solved by the proposed Algorithm jT] Furthermore, it can be shown that 
Pi — p* is a subgradient for Wi, where p* is the solution of (l32l) for the current iteration {m,}. 
Therefore, the solution of (|3TT) can be obtained by the projected subgradient method: for a set of 
given weights j j , we can obtain / ({^ j j an d the corresponding {p *} by the proposed 
Algorithm [I] in the nth iteration; after that, the weights can be updated as 

ru- n+1) = max + t n {P l Vi e S. (33) 

By iteratively updating the weights {w^} and solving (l32l) . we can obtain the solution of (1371) . 
Finally, iterative removal of users with their corresponding SAPs as discussed in Section Ull-B I can 
be applied according to the solution of 071) until all the remaining users can be admitted. 

IV. Asymptotic Analysis and Algorithm Design for Large Systems 

The algorithm proposed in Section [HI] requires instantaneous channel state information. However, 
the instantaneous channel state information may change rapidly with time, which can incur frequent 
resumptions of Algorithm Q] to determine the SAP admission. On the other hand, 5G wireless 
networks are envisioned to be characterized by large numbers of antennas at WBH and dense 
deployments of SAPs. Under those circumstances, certain system parameters tend to become 
deterministic quantities that only depend on large-scale channel statistical information and the QoS 
requirements at SAPs. The large-scale channel statistical information include pathloss, shadowing 
and antenna gain, which do not change rapidly with time. The optimum power allocation p* in 
© and the final selection of SAPs based on x* will also tend to be deterministic irrespective of 
the actual channel changes. In the following, we use random matrix theory to provide an iterative 
algorithm for the SAP admission control problem of large systems. Such an iterative algorithm 
solves © in the asymptotic optimum sense. As long as the large-scale channel coefficients and the 


9iiPi - Vi y ■ ,- <? 

—n -> -, Vi e s 

y..i \.j/i Ui.il'j ■ "i 1 + x i (32) 

E i£S W iPi < Ei G 5 W * P L 

X > 0. 
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QoS requirements remain unchanged, the selected SAPs using large system analysis will almost 
surely be the same as the results obtained by the method of Section HIT] based on instantaneous 
channel state information. 

We assume the number of transmit antennas M at WBH and the number of SAPs N go to infinity 
while the ratio N/M remains bounded, i.e., let M, N —> oo while 0 < lim inf jj < lim sup fi < °°- 
Such an assumption is abbreviated as M —* oo. We use -'"A to denote almost sure convergence, 
where /(x) -^A a means a is the deterministic equivalent of /(x) as M —> oo. The following 
fading channel model is used for large system analysis 

h< = y/dihi, Mi = 1, • • • , N (34) 


where di represents the large-scale fading coefficient between the WBH and the ?’th SAP. The h, 
denotes a normalized channel vector whose elements are independent and identically distributed 
(i.i.d.) CAf(0,1) random variables. 

Under the channel model (]34l) . the uplink SINR ([8]) using MMSE receive beamformers can be 
expressed as 


SINRV (q,{uf MSE ^ 


I) 


M 1 


Qidi{- H 


N 


-1 


y %h,h H 

■4—/ M J ■ 


N 


M~° j 
Qjdj- 


Wil 


+ Wi I 1 hi Mi. 


M h i \ ^ M 


(35) 

(36) 


Calculating the uplink SINR (l36l) involves matrix inversion of dimension M, which becomes 
increasingly complex as M —>• oo. However, the following lemma shows that the uplink SINR (1361) 
will tend to deterministic quantities asymptotically. Such deterministic quantities only depend on 
large-scale fading coefficients {di} irrespective of the instantaneous channel state information { lx,}. 

Lemma 4: As M —)■ oo, the uplink SINRs (1361) using MMSE receive beamformers approach 
deterministic quantities for given q, i.e., 


SINR 


(q, {uf MSE }) 


qidi<pi( q), as M ->• oo 


(37) 


where 


N 


w(q) = + x? E 


Qjdj 


-i 


M j ; 1 + Qjdj<Pi( q) 


Mi = !,-■■ ,N. 


Proof: By applying the trace lemma [f3TTl . If32l and the right-sided correlation model 
the proof follows ll32l Remark 6.1]. 


(38) 


□ 
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The value of </?j(q) is defined implicitly by fixed-point equation (l38l) and it is deterministic for 
given q. 

Similar to the uplink SINR using MMSE receive beamformers, the equivalent channels of the 
system using MMSE beamformers also approach deterministic quantities asymptotically. Using the 
channel model (l34l) . we have 


1 

M 


|hfuf MSE | 2 = ^ 

hfuf MSE 

\ l J 1 M 

1 J 


V*,j = !,•••, N. 


(39) 


Depending on whether i = j, the equivalent channels approach different deterministic values 
asymptotically. The following lemma shows the asymptotic equivalent channels when i — j. 


Lemma 5: As M —> oo, the equivalent channels -f |hfu; 
approach deterministic quantities, i.e., 

1 

M 

where 


MMSE 


h ff u MMSEI 2 a - s -, 


diV 2 M) 

-&( q) 


using MMSE beamformers 


(40) 


N 


v'ii q) = • w i 




qj dj 


-1 


(41) 


j=w (1 + qAvM))' 

Proof: By utilizing the right-sided correlation model ll33ll . 041 and its derivatives, the proof 
follows [[35] Theorem 2], □ 

The asymptotic equivalent channel when i f j can be obtained by the following lemma. 


Lemma 6: As M —)■ oo, the equivalent channels |hfu^ 
proach deterministic quantities, i.e., 


MMSE I 

j 


using MMSE beamformers ap- 


MMSE | 2 a-s-, 


I h >, 


di 


when i f j. 


(42) 


(1 + ^G^-(q)) 2 ’ 

Proof: By utilizing Lemma [5] and applying the matrix inversion lemma P3l and the rank-1 
perturbation lemma |[36l . the proof follows 051 Theorem 2]. □ 

As M —> oo, we can substitute the deterministic equivalents of the system parameters into 
Algorithm Q] and obtain an iterative algorithm for large systems that solves © as shown in 
Algorithm [2] Lemma [4] is applied to (l23l) when updating q. Lemma [5] is applied to (l28l) and 
(l29l) when updating v and p. Lemma [6] is applied to P()l) when updating {z/j}. Algorithm [2] obtains 
the asymptotically optimum solutions for p*, q* and x* of © using MMSE beamformers when 
M —> oo. Note the beamformers {u*} do not approach deterministic quantities because their values 
depend on the actual channel realization. In real system implementations, the outputs of Algorithm[2] 
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Algorithm 2 Large system iterative algorithm to solve © 


l: Initialization: <^;(q) > 0, > 0 and q t > 0, Vi, such that i n iQi — MP. 

2: repeat 

3: Store old values of q 

q = q. 


4: Update 


Qi = 


li 


max(z/j, 1) • </?j(q)dj’ 


Vi = 1 , 


N. 


(43) 


5: Normalize 


MP _ 

q= -^q- 

q 


(44) 


6: Set 


q = ^(q + q)- 


(45) 


7: Update <^j(q) 


^i(q) = 


Wi 


1 


E iV 

3 = 1 )jV ? ' (fi (q) 




Vi. 


Calculate the corresponding </?'(q) 

^i(q) = 


-¥>t(q) 


I 1 

Wi M 


Qjdj 


Vi. 


1 (l+9jrfj¥>i(q)) 


(46) 


(47) 


9: Calculate /i 




—M^7^(q) 
max(z/j, 1) • qiditf(q)' 


(48) 


10 : Calculate the downlink power 


-M^7^'(q) 
max(z/j, 1) • %d^f(q)/C 


Vi. 


(49) 


li: Update z/j 


z/i = n* + 


1 

M 


Pjdi 

(1 + q i d i ip j (q)Y 


Mi 

M ' 


Mi. 


12: until |o,- — oJ < e, Vi. 


(50) 
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can be calculated and stored. Depending on the actual channel state information {hj}, the MMSE 
beamformers {u*} can be directly obtained using ©. Unlike Algorithm[H there is no need to resume 
the algorithm for different channel realizations. In addition, the large-scale fading coefficients of 
SAPs changes slowly in reality. Algorithm [2] can be re-invoked as soon as the QoS requirements 
of SAPs change. Therefore, Algorithm [2] can save lots of computation in large systems. After we 
obtain the output of Algorithm [2] iterative SAP removal can be applied based on the result of x 
achieved by Algorithm [2] 

A. Complexity Analysis 

For the finite system analysis, the iterative Algorithm [j] needs to perform iterative updates of the 
primal and dual variables. The most computational intensive steps in the iterations are the steps 
involving matrix inversion and matrix-vector multiplication in (1231) and (l26l) . where each step has 
the complexity of 0(NM 2 + M 3 ) for each user. Considering all the N users, the complexity of each 
iteration is 0{N 2 M 2 + NM 3 ). By applying the convergence results of fixed-point algorithms [|37l . 
we obtain that the complexity of Algorithm Q] is 0((N 2 M 2 + NM 3 ) log of 1 ), where c\ is a 
constant that determines the convergence speed of iterative operations in Algorithm Q] Following 
similar discussions, we obtain that the computational complexity of the large system iterative 
Algorithm [2] is ()(N 2 log cf ), where c 2 is a constant that determines the convergence speed of 
iterative operations in Algorithm [2] If we use the exhaustive search method, the computational 
complexity will be 0(2 N (N 2 M 2 + NM 3 )N max ), where N max denotes the maximum number of 
iterations using the max-min SINR algorithm ll28l to check the feasibility of the selected SAPs. 
Therefore, the proposed algorithms enjoy polynomial complexity with respect to N, which is much 
lower than the exponential complexity required by the exhaustive search method. 

V. Simulation Results 
A. Convergence Behavior of Finite System Iterative Algorithm 

To verify the iterative algorithm for finite systems, we consider a heterogeneous cellular network, 
where the number of antennas at the WBH is M = 3 and the number of SAPs is N = 4. Each 
SAP has a single receive antenna for the wireless backhaul. The SINR requirement at each SAP 
for the wireless backhaul is set at 7* = 3.01dB, Vi = 1 , • • • , N. The noise variance at each SAP is 
normalized to n* = 1 Watt, Vi. The transmit power constraint for wireless backhaul at the WBH is 


19 



(a) Convergence result of power pi, Vi 



(b) Convergence result of dual variable Vi, Vi 


Fig. 2. Convergence results for finite system iterative algorithms. The solid lines show the convergence of Algorithm |T| and the 
dashed lines show the final results obtained by AUC. 


P = 10 Watt. In the whole simulations, the weights for the transmit power are set to Wi = 1, Vi. A 
randomly generated channel is used to show the convergence result. We compare the results obtained 
by two methods: Algorithm Q] and alternate updates using cvx (AUC). In AUC, we solve © by 
alternatively solving (fTOl) with fixed {uj} using cvx li26l to obtain uplink power q and calculating 
the MMSE receiver {u“ MSE } using the obtained q. In the end, we map the obtained uplink power 
q into the corresponding downlink power p that achieves the same SINR using uplink-downlink 
power mapping Il38l 

The convergence behavior of Algorithm Q] is shown in Fig. [2] Algorithm Q] converges to the same 


result obtained by AUC in less than 8 iterations. Fig. |2(a)| shows the convergence results of {pi}, 
and we initialize all pt with the same value. Fig. |2(b)| shows the convergence results of { v r }. After 
8 iterations, the difference between Algorithm Q] and the final result of AUC is negligible. 

In order to compare the computational complexity of Algorithm Q] and AUC, we perform sim¬ 
ulations on a PC with a Core 2 Duo CPU of internal clock 3.00GHz and 4GB RAM. Twenty 
channel realizations are performed for each algorithm. The number of SAPs varies from N = 5 
to 20. The SINR requirements at the SAPs are set to 7 j = 7.78dB, Vi. The stopping criteria for 
both algorithms are such that the maximum uplink power difference between neighboring iterations 
should be smaller than e = HU 5 . The average simulation time for each channel realization using 
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Fig. 3. Comparison of average simulation time for each channel realization using Algorithm Q] and AUC. 


Algorithm [Hand AUC is shown in Fig. [3j We observe that the average simulation time using AUC is 
5000 to 7000 times longer than that of Algorithm Q] This is because solving the GP problem (fTOl) in 
each iteration requires significant amount of time using standard interior point algorithms compared 
to the iterative Algorithm [H Moreover, the ratio of simulation time increases with the number of 
SAPs N. This shows that the proposed iterative algorithm reduces huge amount of computation to 
solve © without invoking standard convex optimization software. 


B. Convergence Behavior of Large System Iterative Algorithm 

Numerical simulations are carried out to verify the iterative algorithm for large systems. In the 
simulation setup, the WBH has M — 64 antennas and there are N = 32 single antenna SAPs. The 
weights of power are set to Wi = 1, Vi. The power constraint at the WBH is P = 20 Watt and the 
noise variance at each SAP is normalized to n t — 1 Watt, Vi. The SINR requirement at each SAP 
is set to 7 i = 6.02dB, Vi. The large-scale fading coefficients { d t } between the WBH and the SAPs 
are assigned some randomly generated positive values, and they are fixed in the simulation. 

In Fig. gl we compare the results obtained by large system analysis with those obtained by 
Monte Carlo simulations. For the large system analysis, we calculate the transmit power {pi/M} 
and dual values of {z/j} according to Algorithm [2l Those values are the theoretical deterministic 
quantities that {pi/M} and { u ,} converge to when the system dimensions go to infinity, and they 
are shown in vertical bars for each SAP in Fig. |4(a)| and Fig. |4(b)[ respectively. For the Monte 
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dr 3- 
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(a) Transmit powers pi/M, Mi (b) Values of Vi, Mi 

Fig. 4. Comparison of final results in p and v using large system analysis and Monte Carlo simulations. Monte Carlo simulations 
are performed in 100 i.i.d. channel realizations, and their results are represented using error bars. 


Carlo simulations, 100 i.i.d. channel realizations are carried out. For each channel realization, the 
transmit power {p^/M} and dual values of { is, } are obtained using Algorithm Q} For each SAP i, 
the values of Pi/M and //, obtained by Monte Carlo simulations are shown in error bars, where the 
mean value is indicated by a cross sign and the distance above and below the mean value denotes 
the standard deviation. We observe that the mean values accurately match the results obtained 


by large system analysis and the standard deviations are usually very small in both Fig. 4(a) 


and Fig. |4(b)[ This observation shows that {p,/M} and {;/,} really converge to the deterministic 
quantities predicted by large system analysis irrespective of the actual channel realization as the 
system dimensions become large. Instead of the instantaneous channel information, the large system 
analysis Algorithm [2] utili z es the large-scale fading coefficients { d ,} as the channel inputs and 
requires infrequent updates. 


C. Cellular Network Simulations 

We perform numerical simulations for SAP admission control using a macrocell network setup 
with multiple small cells. The WBH is located at the center of the cell. We assume the MBS 
is co-located with the WBH and has the same number of antennas. The cell radius is 1km. The 
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TABLE I 

Summary of simulation parameters 


Number of macrocells 

single 

Cell radius 

1 km 

Small cell radius 

30m 

WBH transmit antenna gain 

5dB 

WBH transmit power constraint 

30dBm 

Log-normal shadowing 

lOdB 

Transmission spectrum for backhaul 

10MHz 

Noise variance 

-93.98dBm 

SAP location distribution 

uniform 

Number of antennas per SAP 

1 

Number of antennas per user 

1 


transmit antenna gain at the WBH is 5dB. The pathloss model from the WBH to the SAPs is 

L(dB) = 128 + 37.6 • log 10 D, (51) 

where D represents the distance between the WBH and the SAP in the unit of km. The log-normal 
shadowing parameter is KMB. The bandwidth of the wireless backhaul is 10MHz. The WBH 
transmit power constraint is 30dBm. The noise variance at each SAP is -93.98dBm. The cell radius 
of each small cell is 30m. We assume the SAPs are randomly and uniformly distributed within the 
cell. The channel pathloss, shadowing parameters, transmit power constraints and antenna gains are 
based on If39| . Table H] shows a summary of the above simulation parameters. 

In the numerical simulations for finite system SAP admission, we compare the proposed SAP 
admission control method, which employs Algorithm Q] and iterative SAP removal, to other two 
commonly used methods in literature: the Lagrange duality based heuristic removal user admission 
(HRUA) algorithm lf40l and the semiorthogonal user selection (SUS) algorithm li4ll . We also 
compare the results with the exhaustive search (ES) method. The ES method obtains the maximum 
number of admitted SAPs by searching through all possible choices of SAPs and chooses the set 
of SAPs with the maximum cardinality. Even though the ES method produces the optimum results 
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(a) Average number of admitted SAPs 


(b) Transmit power yVpi 


Fig. 5. Cellular network simulations for finite system SAP admission. 


for ©, the computational load of the ES method is very high and it is only used as a benchmark 
for comparing different methods. We consider a network where the WBH has M = 4 antennas 
and there are TV = 12 SAPs in the network. We generate 60 channel realizations for each SAP 
location layout and 20 different SAP location layouts are performed in simulations. The results 
for the considered methods are shown in Fig. [5j The average number of admitted SAPs for each 


channel realization is shown in Fig. |5(a)[ We observe that the average number of admitted SAPs 
by our proposed iterative method is nearly identical to that obtained by the optimum ES method. 
The HURA algorithm supports one less SAP compared to the ES method on average. The results 
achieved by the SUS algorithm is the worst, which is 2 to 1.2 less than the ES method. The transmit 


power YhiPi ’ s shown in Fig. |5(b)[ Our proposed algorithm generally has higher power compared 
to HURA and SUS, but all the considered algorithms satisfy the transmit power constraint P. 

We also perform simulations using the cellular network parameters in Table Q] for the case when 
the QoS requirements at different SAPs are unequal and change with time. In the simulation, each 
SAP is responsible for serving one user. We generate 60 different SAP location layouts. For each 
SAP location layout, 20 different user locations in each SAP are considered. In each time slot, the 
SINR requirement of each user is drawn randomly from 4.3dB to 18.7dB. The simulation results 
for different SAP admission schemes are shown in Table [TO After the admission of SAPs, we use 
the max-min SINR algorithm l!T6l to verify whether the users in the admitted small cells can be 
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TABLE II 

Comparison of SAP admission control schemes with unequal SINR requirements 



Ex. Search 

Proposed 

HURA 

SUS 

Average admitted SAPs / users 

3.49 / 3.42 

3.29 / 3.21 

1.78 / 1.75 

1.80 / 1.78 

Average WBH transmit power (dBm) 

28.73 

27.34 

25.66 

24.61 

Average SAP / user sum rate (b/s/Hz) 

10.97 / 10.69 

11.16 / 10.79 

6.92 / 6.80 

6.46 / 6.35 


simultaneously supported for their SINR requirements. If not, we apply the user admission control 
discussed in Section IIII-CI to select the maximum set of users to satisfy their QoS requirements 
within the admitted small cells. The user rates are obtained from the finally admitted users. For the 
ES method, there may be more than one set of SAPs that have the maximum cardinality. In that 
case, we randomly choose one set of SAPs from them due to complexity issues. We observe from 
Table UH that the average number of admitted SAPs for the proposed scheme is very close to the 
optimum results obtained by the ES method, which far outperforms the HURA and SUS methods. 
The average user sum rate achieved by the proposed scheme is also 58.7% and 70% higher than 
the HURA and SUS methods, respectively. In Table H the number of finally admitted users is quite 
close to that of the admitted SAPs. This shows that the small cell user admission control only needs 
to be performed rarely. 

The simulation results comparing the user throughput with and without SAPs in a cellular network 
is shown in Table Hill We use the simulation parameters in Table U and the simulation setup is similar 
to that described for Table HU We compare two cases: one is that the WBH performs SAP admission 
and transmit data to the admitted SAPs using wireless backhaul. Then the admitted SAPs transmit 
those data to their users using in-band channel taking into account the intercell interference. If the 
admitted users cannot be simultaneously supported for the given SINR requirements, user admission 
control discussed in Section IIII-CI is applied. The other is that the users are served directly by the 
MBS using in-band channel. In both cases, the transmit power of in-band channel is the same. 
In Table HUl we observe that the user throughput with SAP admission is about 8.3 times that 
without SAPs. This is because the deployment of SAPs significantly reduces the distance between 
the transmitters and users. Without using small cells, the signal strength at the users will be much 
weaker when transmitted from the MBS. This shows that small cells can boost the system throughput 
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TABLE III 

Comparison of user throughput with and without SAPs 



with SAPs 

without SAPs 

Average admitted (SAPs) / users 

3.31 / 3.21 

0.43 

Average (SAP) / user sum rate (b/s/Hz) 

11.17 / 10.75 

1.30 



< 25 


I 20 


—©— Asymptotic, 30dBm 
x Mean value, 30dBm 
- © - Asymptotic, 27.0dBm 
A Mean value, 27.0dBm 
■ -O- Asymptotic, 24.8dBm 
* Mean value, 24.8dBm 


10 12 14 

SINR requirement (dB) 


-©— SINR requirement 
x Mean value, 30dBm 
A Mean value, 27.0dBm 
* Mean value, 24.8dBm 


£ 12 
CL 


10 12 14 

SINR requirement (dB) 


(a) Average number of admitted SAPs 


(b) Min. SINR for SAPs admitted using large system analysis 


Fig. 6. Comparing results using large system analysis and Monte Carlo simulations. Monte Carlo simulations are performed in 
100 i.i.d. channel realizations, and their results are represented using error bars. 


considerably. 

The simulations for large systems are carried out employing cellular network parameters of 
Table HI In the simulations, the WBH is equipped with M = 64 antennas and there are N = 32 
single antenna SAPs. The weights of power are set to Wi = 1, Vi The power constraints are set 
to P = 30dBm, 27.0dBm, and 24.8dBm. We compare the simulation results using large system 
analysis and those using Monte Carlo simulations in Fig. [6} Fig. |6(a)| compares the number of 
admitted SAPs under given power and SINR constraints. The asymptotic results are shown in 
circles, and they are obtained by the SAP admission control method employing Algorithm [2] and 
iterative SAP removal using large-scale fading coefficients. Those asymptotic results are compared 
to the Monte Carlo simulation results, where 100 i.i.d. channel realizations are carried out. For each 
channel realization, the Monte Carlo simulation performs SAP admission employing Algorithm |Tj 
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and iterative SAP removal. Those Monte Carlo results are shown in error bars. We observe from 


Fig. |6(a)| that the standard deviations of the Monte Carlo results are very small, which shows 
the Monte Carlo results tend to be deterministic values. Such deterministic values are accurately 
predicted by the asymptotic results using large system analysis. In Fig. |6(b)[ we show the achieved 
minimum SINR by the SAPs that are selected employing the large system analysis. For the given 
SINR constraint and large-scale fading coefficients, the set of admitted SAPs is selected using the 
SAP admission control method employing Algorithm [2] and iterative SAP removal. The max-min 
SINR algorithm lf28l is utilized for each channel realization to maximize the minimum SINR for that 
set of admitted SAPs over 100 i.i.d. channel realizations. The error bars show the achieved SINR 
values for different channel realizations, which are very close or above the SINR requirements. 


Fig. |6(b)| shows that the admitted SAPs can satisfy the SINR requirements in nearly all channel 
realizations even though the SAPs are selected using large system analysis. 


VI. Conclusion 

We considered the problem of SAP admission control in a heterogeneous cellular network using 
wireless backhaul. In order to divert the users from the macrocell to small cells, as well as to 
minimize the total cost of building backhaul to serve small cells in the network, the WBH needs 
to simultaneously serve as many SAPs as possible under given power and SINR constraints. Such 
a problem is combinatorial and NP-hard. We applied C-norm relaxation and proposed an iterative 
algorithm to solve the relaxed problem. The local convergence property of the iterative algorithm 
was proved. Based on the solution of the C-relaxed problem, the SAPs were iteratively removed 
until all the remaining SAPs can satisfy the power and SINR constraints. We also proposed a 
large system iterative algorithm using random matrix theory. Such a algorithm only requires large- 
scale channel coefficients to perform SAP admission control for large systems irrespective of 
instantaneous channel information. Simulations showed that the finite system iterative algorithm 
achieved near optimum results and the large system iterative algorithm predicted the Monte Carlo 
simulation results accurately. 
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Appendix A 
Proof of Lemma Q] 

For any set of given beamformers {n,}, according to fl42i Proposition 27.2], both the uplink 
and downlink have the same SINR feasible region under the power constraint P with the uplink 
and downlink SINR definition of © and ©, respectively. The target SINRs 7 ,;/(1 + x t ), Mi, are 
feasible in the downlink if and only if the same targets are feasible in the uplink. Therefore, the 
same set of x is feasible for © and © under the power constraint P. Since this is true for any set 
of beamformers {n,}, the optimum solution of {u,} and x in © and © are the same. For given 
{u,} and x, the mapping between uplink power q and downlink power p that achieves the same 
SINRs can be obtained by the uplink-downlink power mapping [[27]. 

Appendix B 
Proof of Lemma [2] 

Consider the uplink power allocation q*, it is known ll43l that the uplink SINR © is maximized 
by the MMSE receiver, i.e., SINR^ (q*, |u“ MSE }) > SINR^ (q*, {uj}). If u* 7 ^ u mmse j n 
solution of ©, substituting Uj with u“ MSE will also satisfy ©. Therefore, |u“ MSE j must be the 
solution of optimum receive beamformers in ©. The corresponding SINR can be obtained by 
substituting (© into ©. 


Appendix C 
Proof of Lemma [3] 

The Lagrangian of (fTOl) can be expressed as 

7* (V,/, C>."‘ + n\ 


£(q, {Ui} , x) = ^(1-07)07+^ ^ log ■ 


Gael’’ 


log(l + Xi 




According to the KKT conditions, we have 


= 1 - on - V ' = 0, Vi = 1, ■ ■ ■ ,N] 

OXi 1 + Xi 


dC 

dqi 


= ~Vi 




Gije qi 


Yjkrt Gkje qk + Wj 


+ firiie qi = 0, Mi. 


rtie qi — P 

(52) 

(53) 

(54) 
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Therefore, we have 


Xi = max (z/j — 1,0), Vi; 


Vi = 


VjGij 


Jfi 'Ek* j Gk j e*’'+w j 


+ A mi 


oli 


Vi. 


(55) 

(56) 


Because eV > 0, any z/* = 0 requires /i = 0 and z/j = 0, Vj = 1, ■ • • , N, simultaneously. By 
abandoning this trivial solution, we have /i > 0 and z/j > 0, Vi. According to the complementary 
slackness conditions HH and substituting P H = qi/M, we have 


(1 + Xi)Gaqi 

Mli 

If 

i 

Substituting (1571) into (l56l) . we have 




Qj 


M 


+ Wi 


P. 


Vi; 


(57) 

(58) 


Therefore 


We define 


Mi/j 

Qi 


E MGij^jVj 

(TT^- + " n 


Vi. 


(1 ± 3 ±r M i m = sr r MljVj + r> v» 

7i "(1 + xOG^/i ^ '-'(1 -TiV:.uq.il! ' 11 


Pi_ A M^jVj 
M (1 + Xi)Guqin’ 


(59) 


(60) 


( 61 ) 


and we have 

(1 + Xi)GiiPi _ \ - GijPj t _ Wi . 

M 7i _ “ M +n " 


(62) 


Multiply both sides of (l57l) with p, and sum them up for all i. We also multiply both sides of (l62l) 
with qi and sum them up for all i. Because JA Yj^i G jGjPi = Yi Y, y i ( ’i.i , PP.r we have 





niqi _ p 

^ M 


(63) 


Eq. (1621) shows the power allocation p defined in (16TI) achieves the same SINR as the uplink q. 
Eq. (l63l) shows p satisfies the same power constraint as q. When the variables on the right-hand-side 
of (ED) are the optimum primal and dual solutions of ©, the power p defined in (I6TI) corresponds 
to the optimum downlink power that solves © according to Lemma [Q 
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Appendix D 
Proof of Theorem Q] 


We only show the proof assuming that u* > 1, Vi, at the optimum solution of ©. The proof 
that some v* < 1 can be obtained likewise. Since we consider points that are sufficiently close to 
the optimum solution, we assume the beamformers {u,} to be sufficiently close to the optimum 
beamformers {u*}. Then the equivalent channel is Gij = |hfu*| . Vi,j, and it can be assumed to 
be fixed. The updates in Algorithm Q] boil down to the updates of q, and //,, Vi. We define a vector 
u> e and a matrix F e where oji = M r Wi%/Ga, \/i — 1, • • • ,N, and 


-if 1 - if /■ / 

^03 

0 , if i = J. 


(64) 


Furthermore, we introduce a vector y = [yi,--- , yf T , where y t = u,/ q,, Vi. Therefore, u = y o q. 

The updates of q in (l23l) is actually obtained from the KKT condition (fl3l) using MMSE receivers 
{u,}. Since we assume u* > 1 and u, « u* here, the updates of q can be expressed as 


q (m+1 ) = diag (i/ (m) ) 1 (F T q (m) + iv) 

= diag (q (m) o y (m) ) _1 (F r q (m) + w) . 


(65) 

( 66 ) 


By substituting (l28l) and (l29l) into (l30l) . the updates of y, which is obtained by y = v> o q ', can 
be expressed as 


,<m+l> = F . ( q W) 1 + J_ n . (qC”))' 1 ) 

(m)'v — ^ 


= F + 


MP 


-lieu 


q m ) 


(67) 

( 68 ) 


T 

By dropping the time indices and letting z = [q r , y 7 ] , the fixed-point updates of z can be 
expressed as 

diag(q o y) _ 1 (F 7 q + uj) 

( F + ifc n " T ) <r' 


T(z) = 


/i(q, y) 
M q, y) 

Its Jacobian matrix can be written as 

dfi/dq T df 1 /dy T 
df 2 /dq T df 2 /dy T 


(69) 


/ diag(q o y) 1 F i 0 




E- 


0 


0 


diag(qoy) 2 0 

0 diag(qoy )- 2 


(70) 


(71) 
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where 


E = 


diag(F r q + u>) diag(y) diag(F T q + u>) diag(q) 

( F + wp nujT ) dia §(y ) 2 0 

At the optimum solution, we have 

q* = diag(q* o y*) _1 (F T q + u>) 


r=[F+ MP™ )(q * 


Substitute (1751) back to (1721) . Let J* = J(q = q*, y = y*) and A = J* + I, we have 

diag(q* o y*) _1 F T diag (q*) diag (y*)^ 1 
(F + ^naj 7 ) diag (q*)“ 2 I 

The matrix A is nonnegative and irreducible. Furthermore, we have 


A = 


A 


according to (1731) and (1741) . Therefore, 


/qA 

| + ( 

l y *) 

1 \ 


diag(q* o y) 1 uj 

0 


= 2 


(72) 


(73) 

(74) 


(75) 


(76) 


A 


$ 2 


(77) 


and (q*) , (y*) ^ 0. Because A is nonnegative and irreducible, (1771) implies its spectral radius 

P( A) < 2 according to 051 Theorem 1.11]. 

To ensure contraction mapping of the algorithm, the step (1751) must be invoked. Consider the 
update 


(m+l) = I z (m) + Ijv Z M). 

2 2 y ’ 

Let z ^ = z* — e( m ) and z( m+1 ) = z* — e( m+1 \ where z* is the optimum solution. Then 

z* - = i(z* - e (m) ) + —T(z* - e (m) ) 


2 ' 2 

1 '- (m) ) + ^ (T( Z *) - J*-£ (m) ) 


~ - Z — £ 
2 v 


= Z*- ( ^I+ij*)s (m) . 


Here we used z* = T(z*). Therefore, 

(m+l) 




(78) 

(79) 

(80) 
(81) 

(82) 
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Because p( A) < 2, we have p {\I + |J*) < 1. This shows there exists neighborhood of the optimum 
solution, which ensures the mapping in Algorithm Q] to be a contraction mapping that satisfies the 
Lipschitz condition [f46ll . Therefore, Algorithm Q] converges to the optimum solution if the starting 
point is within this neighborhood. 
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